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^ ' Abstract 

^vq \ We propose a non relativistic effective Lagrangian approach to study hadronic atom 

lO ' observables in the framework of QCD (including photons). We apply our formalism to 

^D , derive a general expression for the width of the ■k~^'k~ atom decaying into two neutral 

11/ ' pions. It contains all terms at leading and next-to-leading order in isospin breaking. 

Qv . The result allows one to evaluate the combination qq — a2 of tttt S'-wave scattering 

lengths from vr^vr" lifetime measurements, like the one presently performed by the 

DIRAC experiment at CERN. 



^' The DIRAC collaboration at CERN |T]] aims to measure the lifetime of the tt+tt^ atom 

r-| ! (pionium) in its ground state at the 10% level. This atom decays predominantly into two 

neutral pions, T = T2Tr" + ^2^ + • • • , with r27/r2,r0 ~ 4 ■ 10~^ |^. The measurement of T2ttO 
allows one 0-01 to determine the difference oq— 02 of the strong S'-wave tttt scattering lengths 
with isospin J = 0, 2. One may then confront the predictions for this quantity obtained in 
c^ I standard ChPT [0, ^ with the lifetime measurement, and furthermore analyze the nature of 

spontaneous chiral symmetry breaking in QCD [Q. In order to perform these investigations, 
one needs to know the theoretical expression for the width of pionium with a precision that 
properly matches the accuracy of the lifetime measurement of DIRAC. It is the aim of the 
present talk to derive a general formula for the tt+tt^ atom decay width in the framework 
of QCD (including photons) by use of effective field theory techniques. Our result contains 
all terms at leading and next-to-leading order in the isospin breaking parameters a ~ 1/137 
and (m„ — "m^)^. On the other hand, we expect that the contributions from next-to- next- 
to-leading order are completely negligible, at least for the analysis of DIRAC data, and we 
therefore discard them here. 

In several recent publications |T^, the decay of vr+vr^ atoms has been studied in the 
framework of a non relativistic effective Lagrangian approach - a method originally proposed 
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by Caswell and Lepage [|TT| to investigate bound states in general. This method has proven 
to be far more efficient for the treatment of loosely bound systems - such as the vr+vr" atom 
- than conventional approaches based on relativistic bound-state equations. It allows one 



e.g. to go beyond the local approximation used in |T2|, |T^. On the other hand, we are not 
aware of a systematic investigation of the decay of the 7r"^7r~ atom in this framework. In 
particular, the chiral expansion of the width has not been discussed, and a comparison of 



the corrections found in this framework with the results of |T2|, |T3] has never been provided. 



The main purpose of my talk - based on the papers 0, [Tj] - is to fill this gap. 

The formation of the atom and its subsequent decay into two neutral pions is induced 
by isospin breaking effects in the underlying theory. In the present framework, these are 
the electromagnetic interactions and the mass difference of the up and down quarks. In the 
following, it is useful to count a ~ 1/137 and [m^ — rriu)'^ as small parameters of order 6. 
More than forty years ago, Deser et al. |^ derived the formula for the width of the tt~p 
atom at leading order in isospin symmetry breaking. Later in Refs. |^, |^, this result was 
adapted to the tt+tt^ atoniFl. In particular, it was shown that - again at leading order in 

27r0 

of the difference ao — 02, 



isospin symmetry breaking effects - the width FoS) of pionium is proportional to the square 



rL°o = I aYiao - a,r ; f = {Ml, - Ml, - \Ml,a^fl\ (1) 



At leading order in 5, the momentum p* becomes w 2Mn-+ (M^+ — M^-o) - this is the expression 

used in @, ||. We prefer to use Eq. (|l]), because in this manner, one disentangles the 
kinematical corrections - due to the expansion of the square root - from true dynamical 
ones. In our recent article 0, we derived a general expression for the pionium lifetime, that 
is valid at leading and next-to-leading order in isospin breaking, 

r2.o = ^aV^2(l + ir). (2) 

The quantities A and K are expanded in powers of 5. In particular, it has been shown in 
i thatO 

A = -Tj-ReA+;°° + o(<5), (3) 

62,71 

where Rey4y^7°° is calculated as follows. One evaluates the relativistic scattering amplitude 
for the process 7r"^7r~ -^ tt^tt^ at order S near threshold and removes the (divergent) Coulomb 
phase. The real part of this matrix element develops singularities that behave like |p|"^ and 
In2|p|/M^+ near threshold (p denotes the center of mass momentum of the charged pions). 
The remainder, evaluated at the n~^7i~ threshold p = 0, equals ReA^^^*^". It contains terms 
of order 6^ and 6 and is normalized such that, in the isospin symmetry limit, A = ao — a2- 

^There are a few misprints in Eq. (6) for the pionium decay rate in Ref. B. The correct result is displayed 
inRef. |]. 

■^ We use throughout the Landau symbols 0{x) [o{x)] for quantities that vanish like x [faster than x] 
when X tends to zero. Furthermore, it is understood that this holds modulo logarithmic terms, i.e. we write 
also 0{x) for xlnx. 



Finally, the quantity K starts at order a In a - its explicit expression up to and including 
terms of order 5 is given by 

K = ^^ (ao + 2a2f - ^ (Ina - 1) {2a, + a^) + o(5), A^ = M^^ - M^o. (4) 

In the derivation of Eqs. (0) - (^), the chiral expansion has not been used. It is the aim of 
my talk to show how this result can be derived (details will be provided in a forthcoming 
pubhcation IQ). 

We proceed as follows. First, we display the non relativistic effective Lagrangian for pions, 
as derived from ChPT. Next, we formulate resonance two-channel vrvr scattering theory by 



applying Feshbach's projection technique [0. This method allows one to explicitly reveal the 
pole structure of the scattering matrix element and to obtain the equation for the bound- 
state energy of the tt+tt^ atom. Solving this equation, bound state characteristic of the 
hadronic atoms are given in terms of the couplings in the non relativistic Lagrangian. At 
the final stage, unknown couplings in the non relativistic Lagrangian can be expressed in 
terms of relativistic tttt scattering amplitudes through the matching procedure. Finally, as 
an illustration of our method, we derive Eqs. (Q) - (^. 

The non relativistic effective Lagrangian C = Cq + Cd + ^c + '^s - at the order of accuracy 
we are working here - consists of the free Lagrangian for charged and neutral pions (Cq), the 
"disconnected" piece [Cd) - providing the correct relativistic relation between the energies 
and momenta of the pions - the Coulomb interaction piece {Cc), and the "connected" piece 
(Cs) which contains local four-pion interaction vertices: 

£o = E 4fo-M., + -^)7r„ (5) 

i=±,0 ^ -^^VIt,. J 

^D = Y. ^1(5^ + ---)^- ^c = -47ra(7rL7r_)A-i(7rK) + ---, 

Cs = Ci7r;'^7rl7r+7r_ + C2[7r;'^7rl(7ro)^ + h.c.] + C3(47ro)^ 

+ chilli A 7il{7iof + 7r:^7rl7ro A ttq + h.c] + ■ ■ ■ , 

where u A v = uAv + f An. The coupling constants q are real at 0{a) and are determined 
through matching to the relativistic theory. 

We now formulate the two-channel tttt scattering theory. We denote the full Hamiltonian 
derived from (^ hj H = Hq + He + V, with V = Hd + Hs- The scattering operator T obeys 
the Lippmann-Schwinger equation T{z) = (Hc + V) + {Hc + V)Go{z)T{z). The free and the 
Coulomb Green operators are defined as Go{z) = {z — Hq)'^ and G{z) = {z — Hq — Hc)~^, 
respectively. The pole structure of the T-matrix is predominantly determined by the static 
Coulomb interaction He, whereas V generates a small shift of the pole positions into the 
complex z-plane and will be treated perturbatively. To this end, we use the method developed 



by Feshbach [T^ a long time ago. The T-matrix in our theory describes the transitions 
between charged |P,p)+ = a+(pi)a_(p2)|0) and neutral |P,p)o = cto(Pi)'2o(P2)|0) states, 
where a| denote the creation operators for non relativistic pions. Further, P = pi + p2 and 
p = |(pi — P2) are the CM and relative momenta of pion pairs, respectively. We work in the 



CM system and remove the CM momentum from the matrix elements of any operator i?, 
introducing the notation a(P, q|-R(2)|0, p)b = {2TT)^6^{P){q\rABiz)\p), where A, B = +, 0. 
The operators rAsiz) act in the Hilbert space of vectors |p), where the scalar product is 
defined as the integral over the relative three-momenta of pion pairs. 

In order to avoid the complications associated with charged particles in the final states, we 
consider the elastic scattering process tt^tt^ -^ tt^tt^. In the vicinity of the vr+vr" threshold, 
the scattering matrix element develops a pole at [|T3| 

z-Eo-i^o\r++iz)\^o) = 0, (6) 

where (p|^I/o) = ^o(p) stands for the unperturbed Coulomb ground-state wave function, 
and Eq is the corresponding ground-state energy. According to the conventional definition, 
the decay width is F = —2lm.z. The operator tab{z) denotes the "Coulomb-pole removed" 
transition operator that satisfies the equation 

TAsiz) = Vab + VA+9++iz)T+B{z) + -VAo9ooiz)ToBiz), 

(q|^++(^)|p) = (q|^++(^)|p) — 7 — ^TY- ■ (n 

[Z — _C/oJ 

It remains to solve Eq. (0) in the dimensional regularization scheme and with the use of the 



effective potential technique (see details in p, Q)- We find that the width of the 7r"'"7r 



atom - in terms of the effective couplings q - is given by 



— ^p^/2^^0 1 + — ^ (c2-2pc4)2 l-pc^-^ 1-— ^eci + 



TV ' 

where p = 2M^o(M^+ - M^o - M^+a^/8), ^ = 21na-3 + A + \n{M^+ / /j,^) , 
A = {fi'^Y~^[{d — 3)~^ — r'(l) — ln47r]. The ellipsis denotes higher order terms in isospin 
breaking. The divergent term proportional to A stems from a charged pion loop with one 
Coulomb photon exchange. It is removed by the renormalization procedure in the scattering 
sector. Next, we consider the matching procedure, which relates the effective couplings 
Ci to the TTTT scattering amplitudes for three channels 7r"'"7r~ -^ ir^n^, tt^t[~ —>■ 7r~^7r~ and 
7r°7r° — >■ 7r°7r° evaluated in the relativistic theory ||16|| : 



3M^+ci=4n{2ao + a2)+o{6), 3 M^+ c^ = 27r {ao + 2a2) + o{6), (9) 



3 



A = -—M^+ 



in 



C2C0 , ,2 \ aM\ ( , , M2 



2c2-4A.(c4 + ^M^oj+^^(^l -A-ln -J^)ciC2 



+ o{^) (10) 



Substituting Eqs. (|^), (|T0|) into Eq. (||), we finally arrive at the general formula for the 7r"'"7r~ 
atom decay width in QCD, given by Eqs. (H) - (H). It is our opinion that these equations 
finalize the attempts to calculate the width r27r0 at next-to-leading order, relegating the 
problem to the evaluation of the physical on-mass-shell scattering amplitude for the process 
TT+TT^ -^ 7r°7r° to any desired order in the chiral expansion. Numerical analysis of the tt+tt" 
atom lifetime in ChPT at one loop, including a comparison with recent work in literature, 
was performed recently in Ref. |3]. 



In conclusion, we have evaluated the width r27rO of the tt+tt" atom in its ground state 
at leading and next-to-leading order in isospin breaking. The non relativistic effective La- 
grangian approach of Caswell and Lepage [|lT| appears to be an extremely suitable tool for 



this purpose, that allows one to completely solve this problem. Its usefulness may be seen 
even more clearly for the case of pvr", pK~, rfvr", dK~ atoms, studied in ongoing or planned 
experiments (PSI, KEK, DA$NE), because this approach trivializes the spin-dependent part 
of the problem. 
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